function [A, b] = processing_boundary_condition_2d(A, b, node, bdary, function_g)

%% Neuman

%% Dirichlet
for k = 1:size(bdary,2)
    if bdary(1,k) == 0
        i = bdary(2,k);
        A(i,:) = 0;
        A(i,i) = 1;
        b(i) = function_g(node(i,1), node(i,2));
    end
end

end